On the universality of small scale turbulence 
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The proposed universality of small scale turbulence is in- 
vestigated for a set of measurements in a cryogenic free jet 
with a variation of the Reynolds number (Re) from 8500 to 
10® (max R}, ~ 1200). The traditional analysis of the statis- 
tics of velocity increments by means of structure functions 
or probability density functions is replaced by a new method 
which is based on the theory of stochastic Markovian pro- 
cesses. It gives access to a more complete characterization by 
means of joint probabilities of finding velocity increments at 
several scales. Based on this more precise method our results 
caU in question the concept of universality. 
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The complex behaviour of turbulent fluid motion has 
been the subject of numerous investigations over the last 
60 years and still the problem is not solved |l| . Especially 
the unexpected frequent occurences of high values for 
velocity fluctuations on small scales, known as small scale 
intermittency, remain a challenging subject for further 
investigations. 

Following an idea by Richardson [|j and the theories 
by Kolmogorov and Oboukhov 1^,^, turbulence is usu- 
ally assumed to be universal in the sense that for scales 
r within the inertial range t] ^ r <^ L the statistics of 
the velocity field is independent of the large scale bound- 
ary conditions, the mechanism of energy dissipation and 
the Reynolds number (Re). Here, L denotes the integral 
length scale and rj the dissipation length. 

Besides its physical impacts, the assumed universality 
of the turbulent cascade has gained considerable impor- 
tance for models and numerical methods such as large 
eddy simulations (LES), cf. Finding experimental 
evidence for the validity of the assumed universality is 
therefore of utmost importance. 

The turbulent cascade is usually investigated by means 
of the velocity difference on a certain length scale r, the 
so-called longitudinal velocity increment u(r) 



u(r) = e • [v (x + er, <) — V (x. t)] 



(1) 



where v and e denote the velocity and an unit vec- 
tor with arbitrary direction, respectively. Traditionally, 



the statistics of u{r) is characterized by its moments 
S^{r) = (u"(r)), the so-called structure functions. For 
scales r within the inertial range, the structure func- 
tions are commonly approximated by power laws in r: 



More pronounced scaling behaviour is 



found for the so-called extended selfsimilarity method [|| . 

Experimental investigations carried out in several flow 
configurations at a large variety of Reynolds numbers 
yield strong evidence that the scaling exponents C„ in 
fact show universal behaviour, independent of the exper- 
imental setup 0. A different result, however, was found 
for the probability density functions (pdf) p{u,r). Re- 
cent studies using the theoretical framework of infinitely 
divisible multiplicative cascades show that the relevant 
parameters describing intermittency strongly depend on 
the Reynolds number ||]. 

From the point of view of statistics, a characterization 
of the scale dependent disorder of turbulence by means 
of structure functions or pdfs p{u, r) is incomplete. The- 
oretical studies [|j point out that a complete statistical 
characterization of the turbulent cascade has to take into 
account the joint statistical properties of several incre- 
ments on different length scales. An experimental study 
concerned with the statistical properties of small scale 
turbulence and its possible universalities therefore re- 
quires an analyzing tool which is not based on any as- 
sumption on the underlying physical process and which 
is capable of describing the multiscale statistics of veloc- 
ity increments. Such a tool is given by the mathematical 
framework of Markov processes. Recently, it has been 
shown that this tool allows to derive the stochastic dif- 
ferential equations governing the evolution of the velocity 
increment u in the scale parameter r from experimental 
data |lO|jl|]. 

In this letter we present, firstly, our new method to 
analyse experimental data, secondly, results for different 
Re-numbers, thirdly, experimental findings which ques- 
tion the proposed universality. 

The stochastic process governing the scale dependence 
of the velocity increment is Markovian, if the conditional 
pdf p(ui, ri|u2, 7^2; Mat, rjv) fulfills the relation []l2|,p^: 



p(ui,ri|u2,r2; ...;MAr,rjv) = p(ui, ri|u2, r2). 



(2) 
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The conditional pdf p(iti, ri |u2, r2; uat, r^v) describes 
the probability for finding the increment ui on the small- 
est scale ri provided that the increments U2,...,un are 



1 



given at the larger scales r2,---,rN. We use the conven- 
tions Vi < Ti+i and Ui = u(ri). It could be shown in 
pl| , p^ that experimental data satisfy equation ^ for 
scales ri and differences of scales Ar = r^+i — larger 
than an elementary step size Imar, comparable to the 
mean free path of molecules undergoing a Brownian mo- 
tion. 

As a consequence of (|^), the joint pdf of N increments 
on N different scales simplifies to: 

p{ui,ri;u2,r2; ...■,UN,rN) = r'i|w2, ?'2) x 
xp{u2,r2\u3, r3)...p{uN~i,rN~i\uN, rM)p{uN, tn). (3) 

Equation indicates the importance of the Markovian 
property for the analysis of the turbulent cascade: The 
entire information about the stochastic process, i.e. any 
iV-point or, to be more precise, any A^-scale distribution 
of the velocity increment, is encoded in the conditional 
pdf r|uo, fo) (with r < To). 

Furthermore, it is well known that for Markovian pro- 
cesses the evolution of r|uo, ro) in r is described by 
the Kramers-Moyal-expansion For turbulent data 
it was verified that this expansion stops after the sec- 
ond term pT[ . Thus the conditional pdf p(u,r\uo,rQ) is 
described by the Fokker-Planck equation: 



r-^piu,r\uo,ro) = --^(^D'-^\u,r)p{u,r\uo,ro) 
d 



+ Q^(D^'-\u,r)p{u,r\uo,T^)) . (4) 



By multiplication with p{uo, tq) and integration with re- 
spect to Mo, it can be shown that the single scale pdf 
p{u, r) obeys the same equation. 

Another important feature of the Markov analysis is 
the fact that the coefficients D^^^ and D'-^-' (drift and 
diffusion coefficient, respectively) can be extracted from 
experimental data in a parameter free way by their math- 
ematical definition, see [^ , p^ : 

D^''\u,r)= lim --^MW(M,r, Ar), (5) 
Ar— >o k'.Ar 

M('')(u,r, Ar) = f {u ~ u)^p{u,r - l\r\u,r)du. (6) 



The conditional moments M '^'^^ (u, r, Ar) can easily be 
calculated from experimental data. Approximating the 
limit Ar ^ in equation (ph by linear extrapolation then 
yields estimates for the D^\u,r). 

As a next point, we focus on the analysis of experi- 
mental data measured in a cryogenic axisymmetric he- 
lium gas jet at Reynolds numbers ranging from 8500 to 
757000. Each data set contains 1.6-10^ samples of the lo- 
cal velocity measured in the center of the jet in a vertical 
distance of 40Z3 from the nozzle using a selfmade hotwire 
anemometer [D = 2mm is the diameter of the nozzle). 



Taylor's hypothesis of frozen turbulence was used to con- 
vert time lags into spatial displacements. Following the 
convention chosen in |pT[| , the velocity increments for each 
data set are given in units of ctl = v^cr, where cr is the 
standard deviation of the velocity fluctuations of the re- 
spective data set. 

In order to check consistency of the data with com- 
monly accepted features of fully developed turbulence, we 
calculated the dependence of the Taylor-scale Reynolds 
number R\ on the nozzle-based Reynolds number. Fig- 
ure |l| shows that Rx scales like the square root of Re, 
in accordance with theoretical considerations and earlier 
experimental results. Further details on the experimental 
setup can be found in [|l5|. 
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FIG. 1. Taylor-scale Reynolds number R\ (for details on 
the determination see [11]) as a function of the nozzle-based 
Reynolds number Re. Dotted line: Rx ~ 1.35%/ Re. 

The condition (^) for the Markov property was checked 
using the method proposed in ref. ||T^]. For all the data 
sets, the Markovian property was found to be valid for 
scales ri and differences of scales Ar = r^+i — r^ larger 
than the elementary step size Iman which turned out to 
be of the order of magnitude of the Taylor microscale A 
for all i?e-numbers investigated. 

Having determined the Markov length Imar, the coeffi- 
cients I?*^^^ {u, r) and -D^^^ (m, r) can be estimated from the 
measured conditional moments M'^^^ and M^^^ according 
to equation (|^). The extrapolation towards Ar = was 
performed fitting linear functions to the measured M*^*^^ 
in the intervall Imar ^ Ar < 2lmar [|f6| . 

Figure shows the resulting estimates for the coeffi- 
cients and D^^-' for the data set at Rx = 1180 as a 
function of the velocity increment at several scales r. The 
coefficients exhibit linear and quadratic dependencies on 
the velocity increment, respectively: 



D^^\u,r) = — 7(r)u, 

D''^\u,r) = a{r) - S{r)u + P{r)u^ . 



(7) 



Equation is found to describe the dependence of the 
on u for all scales r as well as for all Reynolds num- 
bers investigated. By fitting the coefRcients D'-'^-' accord- 
ing to (^, it is thus possible to determine the scale de- 
pendence of the coefficients 7, a, i5 and /3. 
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The constant and linear coefRcient of D^^^ a and 5, 
turn out to be linear functions of the scale r (see the inlet 
in fig. I): 



a{r) 



S{r) ^Soj. 



(8) 
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FIG. 2. Coefficients D^^\u,r) (a) and D^^\u,r) (h) as 
functions of the velocity increment u at r = 3X (circles), 
r = L/2 (squares) and r = L (triangles) . The dotted curves 
correspond to linear (a) and polyonmial (b) (degree two) fits 
to the measured data. 

As shown in figure ^, the slopes ao and Sq show strong 
dependencies on the Reynolds number and can be de- 
scribed by power laws in Re with an scaling exponent of 
-3/8: 



ao 



2.8i^e~^/^ So « 0.68i?e-3/^ 



(9) 
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FIG. 3. Coefficients ao (circles) and So (squares) defined 
in eqs. (jj) and as functions of the Reynolds number Re; 
dotted lines represent power laws in Re with an scaling expo- 
nent of —3/8. The inlet displays a{r) as a function of the 
length scale r for Rx = 1180. 



A different result is obtained for the linear term 7(r) of 
Z)'-^\ see figure IJ. It turns out to be a universal function 
of r/A and is found to be well described by 



(10) 



where c = 0.20 ±0.01. 
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FIG. 4. The slope 7(r) of D^^^ (a) and the quadratic coef- 
ficient (3{r) of D'^^^ (h) as functions of the scale r for several 
Reynolds numbers (see legend). 7 is close to a universal func- 
tion of the scale p = r/X (the dotted line is a fit according to 
eq. Aid) ), the coefficient f} exhibits a strong dependence on the 
Reynolds number with a clear tendency towards the limiting 
value f)ao{r) given by equ. ll^ (full line). 



These results allow for a statement on the limiting case 
of infinite Reynolds numbers, Re — *■ 00. According to eq. 
(^) the coefficients a and 5 tend to zero |jl^. 7(r) does 
not depend on Re. Thus drift- and diffusion coefficient 
take the following simple form for Re — *■ 00: 



D^^Hu,r) 
Dg\u,r) 



(3oc{r)u^. 



(11) 



Based on this limiting result we discuss next implications 
for the structure functions S^{r). After the multiplica- 
tion of the corresponding Fokker-Planck equation (|^) for 
p{u, r) with u" from left and successively integrating with 
respect to w, the equation 



nS^{r) 



7(r)-(n-l)/3oo(r) 



(12) 



is obtained. 

According to Kolmogorov's four-fifth law, cf. Q], the 
third order structure function, S'^{r), is proportional to 
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r. Thus, for n = 3, the left side of eq. ( |l^ ) is equal to 
1/3 and f3oo{r) is given by: 



Pooir) 



(13) 



For increasing Reynolds numbers, the experimental re- 
sults for f3{r) in fact show a tendency towards the lim- 
iting value Poo as given by eq. (|l|) (see fig. but it 
is also clearly observed that the convergence is slow and 
that even the highest accessible Reynolds numbers are 
still far from this limiting case. 

To summarize, the mathematical framework of Markov 
processes can succesfully be applied to characterize 
the stochastic behaviour of turbulence with increasing 
Reynolds number. Moreover, the description obtained by 
our method is complete in the sense that the entire infor- 
mation about the stochastic process, i.e. the information 
about any A^-scale pdf p(wi, ri; U2, ri\ ujv, ?'Ar), is en- 
coded in the two coefficients D'^^^{u^r) and D'-^\u,r), 
for which we find rather simple dependencies on their 
arguments w, r and the Reynolds number. 

The i?e~dependence of the coefficients, especially of 
D^'^\ yields strong experimental evidence for a signifi- 
cant change of the stochastic process as the Reynolds 
number increases. This finding clearly contradicts the 
concept of a universal turbulent cascade and might also 
be of importance in large eddy simulations where the in- 
fluence of the subgrid stress on the large scale dynamics 
of a turbulent flow is modeled under the assumption of 
universality. 

It is easily verified that, according to eq. (|lj), the in- 
crease of /3(r) with Re excludes the simple scaling laws 
proposed by Kolmogorov in 1941 j^] even for Re — > oo. 
Furthermore, the universal functional dependence of 7(r) 
on r (eq. (10)) does not support the recently proposed 
constant value of 7 « 1/3 [ |lO| , p^ . The obvious depen- 
dence of the coefficients 7 and /3 on r also contradicts 
the assumption that the structure functions exhibit scal- 
ing behaviour for all orders n, as can be derived from eq. 

With the limiting values for the coefficients D^'''> as 
given by eq. (pi]), the stochastic process for infinite 
Reynolds numbers corresponds to an infinitely divisible 
multiplicativ cascade |jl^ as proposed in ref. How- 
ever, from the slow convergence of the measured coeffi- 
cient (3{r) towards its limiting value (3ooif), it is obvious 
that turbulent data measured in typical laboratory ex- 
periments are still far from that special case. It therefore 
seems questionable to us whether models on turbulence 
established under the assumption of infinite Reynolds 
numbers can be tested in real-life experimental situations 
at all. 
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